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We demonstrate the existence of a large Kondo screening cloud in the k-channel Kondo model 
using both renormalization group improved perturbation theory and the large-k limit. We study 
position (r) dependent spin Green's functions in both static and equal time cases. The equal-time 
Green's function provides a natural definition of the screening cloud profile, in which the large scale 
fif = vf/Tk appears (vf is the Fermi velocity; Tk the Kondo temperature). At large distances 
it consists of both a slowly varying piece and a piece which oscillates at twice the Fermi wave- 
vector, 2kF- This function is calculated at all r in the large-k limit. Static Green's functions 
(Knight shift or susceptibility) consist only of a term oscillating at 2kF, and appear to factorize 
into a function of r times a function of T for rT /vf <^ 1, in agreement with NMR experiments. 
Most of the integrated susceptibility comes from the impurity-impurity part with conduction electron 
contributions suppressed by powers of the bare Kondo coupling. The single-channel and overscreened 
multi-channel cases are rather similar although anomalous power-laws occur in the latter case at 
large r and low T due to irrelevant operator corrections. 
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I. INTRODUCTION 



It is well known that the spin-1/2 impurity interacting antiferromagnetically with a Fermi liquid is completely 
screened at zero temperaturad. This screening is the essence of the Kondo effectH. The question of the screening 
length is much more subtle. Scaling implies, at least dimensionally, that the low energy scale of the model, the Kondo 
temperature Tk ~ Dexp (— I/Aq), should be associated with an exponentially large length scale, = vf/Tk (Here 
Ao is the Kondo coupling times density of states and D is the band width). According to Nozieres' Fermi liquid 
pictures, one could imagine an electron in a region of this size which forms a singlet with the impurity spin. Note 
that this is a more dynamical type of screening than that which occurs for charge impurities in a Fermi liquid since it 
involves a linear combination of states where the impurity spin and the screening electron spin are in either an up-down 
or down-up configuration. In particular, the finiteness of the susceptibility at T ^ should not be attributed to a 
static conduction electron polarization cancelling the impurity spin polarization. Rather it results from the tendency 
of the impurity to form a singlet with the screening electron. 

Whether or not this large screening cloud jieally exists has beeiua controversial subject in the literature, and has 
recently attracted some theoretical interestaQ. Boyce and Slichtertl had performed direct Knight shift measurements 
of the spin-spin correlator at all temperatures and had concluded that there was no evidence of the so-called screening 
cloud. Their measurements, however, were limited to very low distances (not more than several lattice spacings), and 
therefore could not probe directly any possible crossover at the distance scale S,k- 

To study the screening cloud, we will consider the behavior of spatial spin-spin correlation functions, both zero 
frequency and equal time. There are two distance scales in the Kondo problem at finite temperature, and the 
thermal scale = vf/T. On general scaling grounds, the spatial correlators should depend on the ratio of the 
distance r to these two scales. S0rensen and one of uda have suggested a scaling form for the r-dependent Knight 
shift, proportional to the zero frequency spin susceptibility, which has been justified numerically and perturbativelyflfl: 



X(r, T) - p/2 ^ '-p^ f{r/^K, r/^r). (1.1) 



Here we have subtracted the Pauli contribution p/2; p is the density of states per spin. The g-factor for the magnetic 
impurities is not necessarily equal to that of the conduction electrons. This is especially the case for some rare earth 
ions, which have complex multiplet structure. If we take into account this possibility, scaling properties of the local 
spin susceptibility become not sohsimple, and we will consider them below. The Knight shift in this case is a sum of 
two parts, which scale differentlyQ. 

A possible objection to the naive concept of the screening cloud is based on sum rule arguments. The integral of the 



local spin susceptibility Eq.(l.l) is proportional to the zero - frequency correlator {S^iS^^^}, where S^i = J (PrS^i{r) 



1 



is the spin of the conduction electrons, S^^^ = S^i + 5'f,„p. It can be shown that there is no net polarization of the 
conduction electron and this correlator should vanish in the scaling limit (Jp ^ with Tk held fixed). At 
T = this is simply a consequence of the ground state being a singlet. As remarked above, this does not necessarily 
imply the absence of the screening cloud in the sense of Nozieres but only that the screening is a dynamical process. 

In order to see the dynamical cloud of conduction electrons let us consider a snapshot of the system, the equal-time 
correlators. Take K{r,T) = {S^i{r,Q)Sf^p{0)) as an example. Note that (5'f^p(0)S'f„p(0)) = 1/4 for a spin-1/2 
impurity, while (S'|;(0)iS'^q((0)) = as mentioned above. (Note that for this conserved quantity the equal-time and 
zero frequency Green's functions are proportional to each other.) Thus the correlator (S'|;5'f„p) = —1/4; that is 
K{r, T) obeys a sum rule. \K{r, 0)| is a possible definition of the screening cloud profile. 

The ground state properties of spatial correlators are determined by the Kondo scale only. In general we expect 
three different scaling regimes for x(r, T) at a given temperature, with the r-boundaries defined by the thermal and 
Kondo length scales. The goal of this paper is to determine scaling behavior of the spin correlators in these regimes. 

Exponentially large length scale , if present, could have important consequences for the theory of alloys with 
magnetic impurities. Indeed, typical Tk ^ IQK and Ep ~ IQeV makes ^ 10, 000a, where a is the lattice spacing. 
Recently this issue was addressed in ID for Luttingcr liquids with magnetic impuritiesllil, where it was found that a 
crossover happened for riimp ~ 1/^/f. |— ■ 

Although perturbative calculations had beeij-done early onE3, no definite predictions were made regarding the 
size of the Kondo screening cloud. Chen et a/£3 have developed renormalization group approach. They, however, 
only considered short-range correlations r <^ ^k- We use the RG-improved perturbative technique, which cannot 
access lowest temperatures T < Tk- In order to gain some insight into what happens at low temperatures, we also 
consider overscreened Simp — 1/2 multi-channel Kondo effect, where the low-temperature fixed point is accessible 
perturbatively using 1/fc expansion, k being the multiplicity of the bands. A very thorough l/k analysis of the multi- 
channel Kondo effect has been performed earlier by Ganj, who, however, camc_tc|_conclusions opposite from ours. We 
also use the recent conformal field theory approach of one of us and LudwigEII u to calculate the properties of the 
low-temperature, long distance correlation functions and the crossover at ^t- This approach, is valid for all k but only 
for r ^ ^K, T <^ Tk and fails to predict the behavior of the spin correlators inside the screening cloud r £^k- The 
result is nevertheless interesting because, as one could expect, the spin-spin correlators reflect the non-Fermi-liquid 
nature of the overscreened multi-channel fixed point. 

The paper is organized as follows. In Section II we introduce the model and remind the reader how it is transformed 
to an equivalent ID model. Wc also define notations which we plan to use in the rest of the paper and derive the scaling 
equations for the spin susceptibility. Section III provides detailed perturbative analysis of the spin-spin correlation 
functions in the ordinary Kondo model (The Fermi liquid fixed point). Section IV is devoted to the Non- Fermi-liquid 
overscreened large-k case, where it is possible to obtain results for the spin correlators at all temperatures and distances 
using the l/k expansion. We discuss our main conclusions in Section V. In Appendix A we mention a few details of 
our perturbative calculations. Appendix B gives the proof of the vanishing of the uniform part of the susceptibility. 
Appendix C gives results on the overscreened case (/c > 1) at T ^ Tk and r S> ■Jif obtained from conformal field 
theory. Some of these results were presented briefly in Ref. [ ^ . 

II. THE MODEL, RENORMALIZATION GROUP AND SCALING EQUATIONS 

In what follows we consider the standard S'imp = 1/2 Kondo model, 

H = J2 efc4"^ka + JSi,„p • J2 V4"^'/'k'/3, (2.1) 
k k,k' 



and the multi-channel S'imp = 1/2 Kondo model. The Hamiltonian for the Simp = 1/2 k-channel Kondo model also 
includes summation over different channels j: 

H = Y. e'^^k^Vka, + JSimp • ^k"'^^k',73- (2.2) 

kj k.k'j 

Summation over repeated raised and lowpred indiccd is implied. The crucial difference between these two models can 
be seen from the form of the /S-functiond: 

/?(A)=-A2 + ^. (2.3) 
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The flow of the effective coupling is different (Fig. |^) for k = 1 and k The low temperature fixed point of the 
multi-channel Kondo problem is shown to have non-Fermi-liquid naturetJ. At large band multiplicity this nontrivial 
fixed po int becomes acce ssible perturbatively. This difference is not important for the purpose of this section, and we 
use Eq.(2.2) and Eq.(2.3) for both multi-channel and k — 1 models. 

The model is simplified if we assume spherically symmetric Fermi-surface. Indeed, linearizing the spectrum and 
observing that scattering only takes place in the s-wave channel, we can expand the wave functions in spherical 
harmonics: 

V'kaj = -^^/J0aj{k) + (...) 

Ho^ J dfce(fc)Vi"^'(fc)^oa,(fc) + (...) 
H^nt = XoVF J dkdk'4''\k)^i;opj{k') ■ S™p (2.4) 

where e{k) = vp(k — kp) is the linearized spectrum near the Fermi surface, (...) are higher harmonics. Here Aq — pJ 
is the dimensionless coupling constant of the Kondo model , and p = kp/2Tr^vp is the density of states per spin. 
The s-wave operators obey standard one-dimensional anti-commutation relations, 

{4"'''ik),i^0c.,nik')} = S:iSf^S{k - k'). (2.5) 

We define left and right movers on a band of width 2 A around kp: 

.A 

i^LM^ / dke^^'^^'Mk + kp). (2.6) 

J -A 

The 3D fermion operators are then written in the form: 

vl.(r) = [e-^'-'-Mr) " e^'^^i^n{r)] + (...), (2.7) 

2v27rr 

where (...) are higher harmonics. The left and right-moving fields defined on r > obey the boundary condition: 

^l(O) =^i?(0). (2.8) 

Flipping the right-moving field to the negative axis, V'l(~^) = V-'-r(^)5 we rewrite the ID Hamiltonian in terms of the 
left-moving field only: 

/CO ^ 
driPl{r){id/dr)4,L{r) + 2^«;^Ao?A1(0)-^l(0) • Si„,p. (2.9) 
-oo ^ 

The purpose of this paper is to analyse various spin-spin correlation functions. The most important of them is the 
distance-dependent Knight shift, which can be measured in NMR experiments. If the impurity spin has a different 
gyromagnetic ratio from that of the conduction electrons, the uniform magnetic field couples to the spin operator 
Sh = Sei + {gs/2)Simp, where Si,„p and Sei = (1/2) / dr1p^r)a■tp{r) is the total spin operator of the impurity and 
conduction electrons, defined with channel sum for the multi-channel problem. The expression for the Knight shift 
then consists of the electron and impurity contributions: 

X(r) ^ [d^< Sli{r,T)Sm Xei{r) + |x.™p(0- (2-10) 

(We set pb — 1-) We will also consider the equal-time spin correlator K{y), defined by: 

K{v) EE< Sl,{r,T^ Q)St„,,{Q) > ■ (2.11) 

The above ID formalism allows to simplify this expression for large rkp ^ 1. Substituting Eq.(2.7) in Eq.( 2.10| ), we 
get: 



^A{r,T)^—^ cos{2kpr) + —^ — , (2.12) 

An r'^vp on r^vp 
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where A corresponds to imp or el. For K(r, T) we get a similar expression: 



i^(r,T) = ^f#cos(2fc,.) + .^«"^^) 



Vf 



(2.13) 



The total electron spin in ID is: 



Scl — 1^ — 

2tt 



dri}l{r)^iJL{r). 



(2.14) 



The uniform and 2kF parts take the form: 



/5 



XunAr,T) = VF j dT < [^lir,T)^Mr,T) +^l{~r,T)^M-r,T)]Sm > 
X2k,,A{r,T) = -VF dr< ^.[(r, r)y ^L(-r, r)51(0) > . 



(2.15) 



Expressions for Kun and Kimp are analogous to those for Xun,imp and X'2kp,imp in Eg. ( 2. 15 ), although they don't 
involve integration over r. 

If the spins of the impurity and conduction electron have equal gyromagnetic ratio (gs = 2), the operat or Si is the 
total spin of conduction electrons and impurity, and is conserved. The Knight shift is then given by Eg. ( 2. 15 ), with 
A = tot. Since the Kondo interaction is local, only boundary (r = 0) operators have non-zero anomalous dimensions. 
Thus the conduction electron spin operator Sei{r) also has zero anomalous dimension, for r ^ 0. The local spin 
susceptibility then obeys the following RG equation: 



D 



d_ 
dD 



-/3(A) 



d_ 

dX 



x{T,X,D,rT/vF) = 0, 



(2.16) 



where D is the ultca- violet cut-off (the bandwidth), and /3(A) is the /3- function. S0rensen and one of us have recently 
made a conjectionti, supported by perturbative and numerical results, that in the scaling limit, rkF S> 1, T ^ Ef, 



■ uitca-v: 
tionEl, sr 

the spin susceptibility has the following form: 



VT 
X 1 — , 

VF 



T 



X2kf 



( rX. JL\ 

\ VF ' Tk J 



47r2 



r^VF 



cos(2fci?r) 



(2.17) 



where X2feF ^ universal functions of two scaling variablesEj. This form follows directly from Eqs.( 2T^ ,2.16). In 
general, one expects that there could be a non-zero phase in Eq.( 2.17] ), and a uniform term. It is easy to see that 
the phase is zero due to particle-hole symmetry^!. Indeed, under particle- hole transformation ipLi^) a'^ip]^{r), so 
Stot S(ot, ip\^{r)(T^pL{—r) — > ip\^{—r)(T^pL{r). Particle-hole symmetry of Eq.( 2.17 ) then requires that the phase 
is zero. This is not so for more realistic Hamilto nian s, for which the particle-hole symmetry is broken. For such 
Hamiltonians there is an additional phase (f> in Eq.(2.17), but this phase does not renormalize. That is, it is essentially 



constant in the scaling region (fc^rr 3> 1). The fact that the uniform part of the spin susceptibility is zero is less trivial. 
For the static local spin susceptiblity we have provedQ that all graphs in perturbation theory contain certain int egral s 
that vanish. These properties hold for the electron and impurity parts of the local spin susceptibility Eq.(2.1f;) 
separately, for both single-channel and multi-channel Kondo effects (see Appendix B). The uniform part and the 
phase are zero for the Knight shift in case of nontrivial gyromagnetic ratio for the impurity s pin (g g 7^ 1) as well. 

Since we consider the problem perturbatively, it is useful to express the scaling function Eq(2.17) in terms of some 
effective coupling constants at an energy scale E, Xe- This way we eliminate non- universal Tk- The energy scales of 
interest are the temperature T and the distance energy scale vf/t. We will denote corresponding effective couplings 
as At and A,. . Expressions in terms of effective couplings can be easily converted into those in terms of Tk , and vice 
versa, provided that the /3-function is known up to the order needed. Indeed, 



dXf 



dlnf 



(2.18) 



where D — wf/A is the bandwidth. Therefore, for the effective coupling at two different energy scales E and E' we 
have: 
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A' 



dX _ E' 



Since Xt^ = 1 can well be regarded as one of possible definitions of Tk , we have 



dX 
W) 



E 



(2.19) 



(2.20) 



and the arguments of the scaling function in Eg. ( ^.17 ) can be replaced by corresponding effecti ve co uplings. 

The renormalization group equations for various parts of the local spin susceptibility in Eq.(2.15) are less trivial. 
Consider first Ximp- Since the Kondo interaction is at the origin,_the fermion bilinear operator has zero anomalous 
dimension, while the operator Simp receives anomalous dimension,LiJ "fimp — A^/2. Rcnormalizability implies that the 
functions XB.imp (B — 2kp, un) obey equations of the form: 



^^+/3(A)^+7™p(A) 



XB,imp{T, X, D, rT/vp) = 0, 



(2.21) 



where 7imp(A) is the anomalous dimension, which in this case is equal to the anomalous dimension of the impurity spin 
operator. The other correlator, XB.eh contains the total conduction electron spin operator Sei, in which integration 
over the electron spin includes a potentially dangerous region near the impurity site. In this region operator mixing 
occurs between the electron spin and impurity spin. Thus XB,ei obeys non-trivial mixed RG equation: 



d d 



XB.ei{T,X,D,rT/vF) = ltmp{X)xBSmp{T,X,D,rT /vp)- 



(2.22) 



This equation can be obtained by subtracting Eq.(2.21) from Eq(2.16) 



shift for (7s 7^ 1 hi terms of Ximp and Xtot, which obey ordinary scaling equations Eqs.( 2.21 
The solution of the scaling equation for XB,imp Eq.( ^.21 ) has the following form: 



It is more con venie n t to express the Knight 

(ill)- 



XB,imp \ Ao, — 



T rT 
An, t:^: 



K Vp 



Z^'o ^<^) nB,^mpiXT,—j 



(2.23) 



^B,irnp { -^T, — \ e 



Jo 



-d\ 



Here b .impiXr , fT / v p) ^ ^B,imp{XTTfT/vp) are some scaling functions to be determined below; Aq = pJ is the bare 
coupling constant. The solution of the scaling equations for Ximp is a function of T/Tk and rT/vp, up to some 

non-universal coefficient. We see that the non-universal coefficient expl— (iA(7imp(A)//3(A))] is equal to unity in 
the scaling limit of zero bare coupling Aq 0, if ^imp{Xo) / I3{Xq) is non-singular in this limit. This is indeed the case 
for the K ondo model. The scaling function <^B,imp{T /Tk^tT jvp), of course, can differ fr om "x ^kwjT /TK.rT /vp) in 
Eq.( ^.17 ). The equal-time correlation functions also obey analogous scaling equations Eq.( 2.21 ) with the anomalous 
dimension which is a sum of the dimension of the corresponding operators (for Ximp it is again ^imp) ■ These equations 
are also applicable to the uniform part of the correlator, which is now non-zero. 

In the rest of this paper we will consider these scaling functions in various regimes, which we now outline. Scaling 
form is applicable for r ^ 1/kp, a nd T ^ _D w Ep. For the single-channel Kondo model perturbative treatment is 
only valid for T ^ Tk- From Eq.(l.l) one could expect that there could be two crossovers: one at and one at 
^K ^ £,T = vp/T. The latter crossover, however, does noi happen as a function of r for T ^ Tk- The lowjtemperature 
correlation functions in the single-channel Kondo effect can be studied using the Fermi liquid approachH. The region 
of validity for this approach is r 3> ^x, T <ti Tk- It provides important information about the low-temperature long- 
distance form of the correlation functions, and the crossover at S> £,k, but is unable to access the most interesting 
region r ^ £^k, and answer the question of existence of the screening cloud. For the multi-channel Kondo effecL 
low-temperature long distance correlation functions can be obtained using the conformal field theory approachllZI' 
which is a generalizarion of Nozieres' Fermi liquid picture. It is also limited to r 3> S,k, T <^ Tk- The interesting 
low temperature region with r ~ ^k only becomes accessible at large k, when the whole scaling function can be 
constructed. 



III. THE SINGLE-CHANNEL KONDO MODEL. 
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A. The Knight shift. 



In what follows we consider the Knight shift in the single channel Kondo model. As mentioned above, the loca|l 
spin susceptibility only has the oscillating part. We have calculated it up to the third order in perturbation theoryQ. 
Summing all the relevant diagrams (see Appendix A), we obtain: 

X2k, [x^—A^.D] = [Ao + XliHD/T) + M{x) + x) (3.1) 

\ Vp J 4smh(27ra;) 

+ Ag(ln2 (D/T) + \n{D/T){2M{x) + 2a; - 0.5) 

+ {M{x) +x){M{x) +0.5) + const)], 



where 



M{x) = In [1 - exp(-47rx)] . (3.2) 



Substituting this expression in Eq.( 2.21 ) we find that scaling is indeed obeyed. At small r, a; <C 1, Eq.( |3.l| ) is rewritten 



X2k. {r, Xa,D)^^ [Ao + A^ In(Ar) + Ag In^ (Ar) + 0.5X1 In(Ar) (3.3) 
- Xl \n{D/T) + Xlconst], 



where A = A-kD/vf = 47rA ^ kp. It is clear from Eq.( |3.3| ) that the infrared divergences of the perturbation theory 
are not cut off at low T by going to small r, as was first noticed by Gancl. In the third order, these divergences 
are associated with the graph shown in Fig. 2. Due to the non-conservation of momentum by the Kondo interaction, 
the bubble on the right gives a logarithmic T-dependent factor which is independent of r. Thus, the interior of the 
screening cloud does not exhibit weak coupling behavior. 

It is convenient to rewrite this result in terms of effective couplings at the energy scales T and vp/r. One can easily 
write down the effective coupling constantE3 at some energy scale oj using the well-known /3-function Eq.(p.3[): 



A^ = Ao + A^ ln{D/uj) + Xl[\n^{D/uj) - (1/2) \n{D/uj) + constant]. (3.4) 

We find that the expression for X2kF is simplified when we use effective couplings At and As at the energy scales T 
and E{x) = T/[l — exp(— 47rx)] e, x = r/^T- When r <C the latter becomes the effective coupling at the distance 
scale r, since E{x) oc vp/r. Eq.( |3.lD in terms of these effective couplings takes the form: 

rT ^ \ {XE + {i'^/'2)X\x + constX\){l~ Xt) 

— I — — 

vp^ J (4/7r^) sinh(27ra;) 



X2kp \X ^ —,Xt ] = / 2A„:-U/o X ' (3-5) 



It is instructive to consider various limiting cases for the scaling function Eq.( |3.5D . For r <C we find: 

X2kAx,XT) = /Sx){Xr + constXl){l - Xt). (3.6) 
If r 3> ^T, the spin susceptibility takes the following form: 

X2k, (x, At) = (3^V4)A|(1 - AT)e"'"^ (3.7) 
For high temperatures T ^ T k there is no crossover 'At r in the behavior of the local spin susceptibility. The 



factor (1 — Xt)/AT in Eq.(3^) is, to the order under consideration, precisely the total impurity susceptibility, xtt {T) . 
This is the total susceptibility less the bulk Pauli term and its value has been determined accuratelytl. Thus Eq. (p!q) 
can be written: 

, , (Ar + constXl) 
X2kAT,r)^ — — ^ — Xtt r , (3.8 

We can compare this result with the experiment of Boyce and Slichter§, who have measured the Knight shift from 
Cu nuclei near the doped Fe impurities, at distances up to 5-th nearest neighbor. At these very small distances of 
order of a few lattice spacings, they have found empirically that the Knight shift obeyed a factorized form, x(r, T) « 
/(r) / (T + Tk), with rapidly oscillating function /(r) for a wide range of T extending from well above to well below 
the Kondo temperature. Although our condition r 3> 1/kp is not satisfied in this experiment, this form coincides with 
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Eq. (3.. 



since the Bethe Ansatz solution for XttiT) niay be quite well approximatecfl by 1/ (T + Tk) at intermediate 
temperatures T ~ Tr-. As one can see from Eq.( p.5| ), this factorization breaks down at r ~ vp/T. 

At low temperatures T <^ Tk and large distances r ^ S^k the behavior of x(r) is determined by the zero-energy 
Fermi liquid fixed pointd. The Kondo impurity acts as a potential scatterer with a phase shift 7r/2 at the Fermi 
surfaced. The local susceptibility follows directly from the formula for Friedel's oscillation in the electron density for 
an s-wave scatterer and 7r/2 phase shift, 



n{r) — no 



1 



■ cos[2kFr + 7r/2] 



27,2^3 

Since the magnetic field H simply shifts the chemical potential by ±g^BH/2 for spin up or spin down electrons, 



(3.9) 



1 dn 
vp dkp 



P 



1 



:Cos{2kpr). 



This implies for the scaling function Eq.(2.21 



X2kf 



The finite-temperature properties of X2kF{r), and, in particular, the croiSjsaver at r 
from the Nozieres' low-energy Hamiltonian for the Fermi liquid fixed pointLrllZI; 



^0 = dri:l{r)^Mr)+^-^Sl,ir), 



(3.10) 
(3.11) 

can be obtained directly 
(3.12) 



where Sei{r) = ip]^{r){(T /2)ipL{r). This definition of Tr- differs from one in Eq.( 2.2C ) or xtt oc 1/{T + Tk) by numerical 
factors 0(1) (Wilson ratio). The expression for X'ikpi'''), Eq.(3.11), is zero order in the leading irrelevant coupling 
constant 1/Tk, and the finite-temperature form of X2kF{r) is easily obtained: 



X2kF{x) 



2ttx 



sinh(27rx) 



rT 
Vp ■ 



(3.13) 



We can derive corrections to Eq. (|3.13 ) by doing perturbation theory in the leading irrelevant operator. For the first 
correction we obtain: 



SX2kpix) 



Tk sinh(27rx) 



(3.14) 



The first correction does not alter the leading order behavior. At zero temperature the scaling function for r ^ 
takes the following form: 



X2kF{r/£.K) = 1 + 



2r ■ 



(3.15) 



This correction gives rise to the first term in the large-distance expansion of our scaling function X2kFij' I ^k)- The 
behavior of the scaling function X2kF in different regimes is summarized in Fig.(^). X2kF{^/^K,r/£^T) exhibits a 
crossover at low T, when the "screening" cloud is formed. At high temperatures this crossover is absent. 

What happens when the (/-factor of the impurity is anomalous? X2kF {f, T) is a sum of impurity and electron 
parts, X2fcF.im p and X2kF.c\- As we have discussed in the previous section, the latter obeys a complicated mixed RG 
equation, Eq.( 2.2^ ). It is more convenient to express the spin susceptibility in terms of the correlators X2fcF,tot(^, 7") 
and X2fcF,imp('', 2^), for which RG equations are simple: 



X2kF {r, T) = (g5/2)x2feF,imp(?', T) + X2kF,c\{'r, T) 

X2kF,c\ — X2kFXot{r,T) — X2kFAmpir,T) 



(3.16) 



Since we have already determined X2/cF,tot('', T), it is sufficient to consider only X2kF-imp{f^ T). From the perturbative 
analysis (see Appendix A) we obtain: 



X2k 



F ,imp 



rT 

X = — ,\,D] = 

Vp 



4sinh(27ra;) 



[Ao + Ag(ln(D/r) + M{x) + 0.5) + 



(3.17) 



Xl{\n^{D/T) + 2\n{D/T)M{x) + {M{x) + Q.bf + const)], 
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where M{x) is the same as in Eq.(3.2). One can easily check that Eq.( ^.21 ) is obeyed with@fl: 



/3(A) = -A2 + ^, ^i^p(A) = y. (3.18) 



We then obtain for the non-universal factor in Eq.(2.23): 



rio -Ti^^^ ^1 + ^, (3.19) 



and the local impurity spin susceptibility takes the following form: 

X2fci.,imp - (^^ + y) ^2fel(AT,a;), 

where the scaling function 



(1)/,, N {\e + constXl){l ~ Xt) , 

X2kA^T,x)^ . . — T (3.21) 

^ (4/7r^) smh(27ra::) 



(3.20) 



differs from that for the conserved local susceptibility Eq.( |3.5| ). For the electron part we obtain: 

^ Ao (1),. , , (37r/2)A|a:(l - At) 
X2..,ei - -yX,..(AT,x) + ■ (3.22) 

The second contribution does not vanish in the scaling limit Aq 0. However, it only becomes substantial at large 
distances r ~ ^t, where there is no additional smallness associated with the factor x = r/^x ■ We conclude that two 
different scaling functions are present in the experimentally measured Knight shift, and their share depends upon the 
gyromagnetic ratios for the impurity and the conduction electrons. 



B. Integrated susceptibilities 

It is instructive to consider the integral of x(r,T) over all space. This quantity determines the polarization of the 
screening cloud in external magnetic field. We immediately see that the contribution from large distances vanishes 
because of the oscillatory behavior of x{r, T) at large r. Nevertheless the integral can be finite due to the contributions 
at small distances r ~ 1/kp. We will specify three different spin correlators: 

Xtt[J- ) = j; , XtiiJ j = 7^ , 

Xii(T) EE / < 5f^p(r)5f^p(0) > dr. (3.23) 



For this choice of correlators the RG equations are simplified, and have the form Eq.(2.21). It seems more natural to 
define correlators of the impurity spin and the total conduction electron spin S^i instead of Stot- 

Xec(T) = r < ^eM^)^cl(O) >dT-xo (3.24) 







Xci(r)^ r<5,^i(r)5f„p(0)>dr, 
Jo 

where xo is the free electron susceptibility, proportional to the volume of the system. However, for this set of spin 
correlators the RG equations are mixed. 

Two of the three spin correlation functions can be measured. The first one is the bulk susceptibility, Xtt{T). The 
electron spin polarization in the p resence of an impurity is determined by the spatial integral of x(^) measured in the 
Knight shift experiment Eq.( 2.1C| ), or, equivalently, by Xtt{T) — Xti(7^)- If the gyromagnetic ratio for the impurity is 



different from 2, the experimentally measured magnetic susceptibility is (s'f/4)xii + QsXie + Xeej while the integrated 
electron susceptibility is given by {gs/2)xie + Xe 
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Since S^^^ is conserved, the spin susceptibilities obey the RG equation, Eq. (2.21), with anomalous dimensions 



determined by the dimension 7imp(A) of the operator S'j^^p. For the three different susceptibilities: 7tt 
and 7ii = 27iiiip. The solutions of these equations take the form Eq.( 2.23 ), 



0, 7ti = 7ii 



(3.25) 



where j labels tt, ti or ii. From our third-order perturbativc analysis using Wilson's resultO for Xtt{T) we have 
obtained that the functions $j(AT) ~ 1 — At coincide for all three susceptibilities up t o and including terms of 
order If this is indeed the case in the Kondo model, we then obtain from Eq.( 3.25|) that in the scaling limit 
Ao ^ both XeeiT) and Xie{T) vanish. At fin ite b are couplings these susceptibilities also become finite, with non- 
universal amplitudes. We then obtain from Eq. ( 3.25| ) for the impurity-electron and electron-electron pieces of the spin 
susceptibility: 



-YXtt(T) 
|xtt(T), 



(3.26) 



Thus, the integrated distance-dependent Knight shift obeys: 

x{r,T)dr = Xee{T)+x^e{T) 



Ao 



Xtt(T). 



(3.27) 



The major contribution to Eq.(3.27) comes from the electron-impurity correlator. It should be emphasized that the 
result is non-zero at finite bare coupling Ao. (A typical experimental value of Ao might be 1/ \ti{Ep /Tk) ~ .15.) It 
is easy to see that the integral in Eq.( |3.27 ) is dominated by r ~ l/kp- Thus most of the small net polarization of 
the electrons in a magnetic field (with the free electron value subtracted) comes from very short distances. However, 
this should not be interpreted as meaning that the screening cloud is small as can be clearly seen from the equal-time 

correlation function discussed in the next sub-section. 

If the equality of the scaling functions ^^(At) defined in Eq. ( 3.251 ) holds at all T, the integrated electron spin 
susceptibility vanishes in the scaling limit of zero bare coupling at all T . The fact that Xie ^^nd Xee are suppressed in the 
scaling limit has been known or conjectured from a variety of different approaches over the years. The earliest result 
of this sort that we are aware of, in the context of the Anderson pipdel, predates the discovery of the Kondo effect 
and is referred to as the Anderson-Clogston compensation thcorcmE£l. It was later established at T=0 from the Bethe 
ansatz solutionEJ. A very simple and general prooill3 of this result follows from the abe lian bosonization approached. 
Beginning with left-moving relativistic fermions on the entire real line, as in Eq. (|2.9| ), we may bosonize to obtain 
left-moving spin and charge bosons. The charge boson decouples and the Hamiltonian for the Kondo Hamiltonian 
can be written in terms of the left-moving spin boson, (f>L which obeys the canonical commutation relation: 



The Hamiltonian becomes: 



H ^ dr 

J —OO 

= 27rvFXQ 



Xr'),drMr)]^(^/mr-r'). 



V 27r 

^l^Ip^ + const ■ (5+e^^^-(°) + h.c.) 
V 27r 



(3.28) 



(3.29) 



Here hi and /ig are the magnetic fields acting on the impurity and the conduction electrons, correspondingly. These 
fields may differ by the ratio of corresponding g- factors. We can get rid of the j drdripLir) term by shifting the 
bosonic field: 

ipL{r) = (pL{r) + hev/vF^^ (3.30) 
The Hamiltonian in terms of the new bosons takes the form: 

hlL 



H 



drvF{dr(pL{r)) 



Attvf 



2TrvFXoS^dr'PLiO) 

heXo 



(3.31) 



+const ■ 2^Ao(^+e'^^^("' + h.c.) - (^h, - 



5^ 
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Thus, our original Hamiltonian with non-zero field he acting on the conduction electrons is exactly equivalent to the 
one with no field acting on conduction electrons and modified impurity field. The same argument was given in Ref.[ 
p5| except that the field shift by heXo/2 was not obtained because another, non-commuting, canoncial transformation 
was performed first to eliminate the z-component of the Kondo interaction. 

In terms of the free energy, this is written as; 

F{h,,he) = -^+F{0,h,~^) (3.32) 
Taking magnetic field derivatives, we easily find: 

'dh 



= -IJIT (3-33) 



d^F _ _Ao 



9/i2 V 2 



where xo = 2mi Pauli term. It is easy to see that this is valid for the anisotropic Kondo model as well, with 

Ao being the z-component of the Kondo interaction. 



C. Equal-time spin-spin correlator 



The equal-time spin correlators provide a snapshot of the Kondo system. The quantity of interest is: 

i^(r,T) = (^|,(r,0)^f„p(0)). (3.34) 
As we have shown in the previous sections, it satisfies a non-zero sum rule: 



/ 



dTK{r,T) = -1/4. 



(3.35) 



The proof that Xun (^) = is based on the fact that the time integral for the Feynman diagrams is zero in all orders 
in perturbation theory (see Appendix B). For the equal-time correlator we don't integrate over the time variable, so 
the uniform p art does not have to vanish. K{r, T) can be rewritten in ID in terms of the uniform and 2kF parts, 
Eqs( 2.17 , 2.15| ). For t he sa me reason as for the impurity p art o f the Knight shift, the equal-time correlator obeys the 
scaling equation Eq.( ^.21 ), with solutions of the form Eq.( 2.23| ). Since the decomposition of K{r,T) into the uniform 
and 2kp parts is only valid in the scaling region kpT 3> 1, the sum rule Eq.( |3.35" ) do cs no t necessarily extend to 
Kuni'f, T). The region r ~ could produce a large contribution to the sum rule Eq.(3.35). 

Consider now the equal-time correlators Kun (t, T) and K2kF T) perturbatively. In the third order we obtain (see 
Appendix A): 



^\-Xl + {l/2)\l-2\lHD/T])T 



exp(27rr/^T) — 1 
r f /, ^ ln(l-e-2-/«T) 



(3.36) 



4sinh(27rr/^T) 



Gi^2{x) are some functions which can be represented as integrals: 



Gi(x) = 
G2{x) = 



2ds 



ds- 



1 1 



In 



(1 - s) (1 - se-4^^) 



■In 



^1 

1 — se 



se 

2irx 



27TX 



(3.37) 



1 - e- 



It is easy to check that Eq.(2.21) is satisfied for both uniform and 2kp parts. The solutions are found in the form 
Eq.( ^.17 ) with the non- universal factor Eq.( [3.19| ). The scaling functions are easily obtained from Eq.( ^.3(: ). The final 
expressions are simplified in the most interesting limiting cases. For r <C ■^t we obtain: 
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Kun{K,r/iT, Ao) = (3.38) 

Zr 

„ /c \ ^ 7r?;i7'Ar-(l + Ao/2) 

8r 

In case of r ^ ^t, these functions take the form: 

K^n{\r.rliT, Ao) = -vr^TAl [\ + ^ g-^-ZC^ (3.39) 

K,,AK.rlir. Ao) = ^ (1 + y) (1 - \Ty-'^'l'^. 

Note that 1^2*;^ is suppressed in this Umit by the small value of the bare coupling. Like the local spin susceptibility, 
the equal-time correlator does not have crossover at r ~ at high temperatures. Instead, the corresponding scaling 
function for r > has a factorized form, K{r / £,t,T /Tk) oc h{r /iT)f2{T /Tk). 

The behavior of the equal-time correlation function at T ^ Tk and r 3> can be calculated using Nozieres Fermi 
liquid approach. Indeed, the impurity spin at the infrared Kondo fixed point should be replaced by the local spin 
density J(0) for r « 0, up to a constant multiplicative factoiOE3: 

t^f Jl(0) .„^ 
Simp oc — . (3.40) 

Substituting this in the definition of K^n and K2kp, we obtain at finite T: 

K2kAr/(T) = -{l/2)KUr/(T) = ^ ■ ^2, ^ , , (3-41) 

1 K smh (ttj r /vp) 

Thus, at T ^ the equal-time correlator K decays as sin kpr/r^ (see Eq.(|]l|)). This result was obtained by Ishii0 
in the context of the Anderson model. The behavior of the equal-time correlators K2kFi'r) and -R',m(r) in different 
regimes is summarized in Fig.^ and Fig.^. 



IV. LARGE K MULTI-CHANNEL KONDO MODEL. 



The information that one gets for the single-channel Kondo model using perturbative RG is very limited, and further 
numerical analysis is required. To justify the presence of the Kondo length scale more, we analyse the multi-channel 
model with large band multiplicity. The generalization of the above perturbative analysis to the multi -cha nnel case 
is quite straightforward. The Hamiltonian for the S'imp = 1/2 k-channel Kondo model is given by Eq.( ^.2[ ). Further 
analysis of Section II applies to the multi-channel jcase as well. Some of the relevant perturbative 1/k calculations for 
the multi-channel Kondo effect were done by Garu. His scaling equations and conclusions about the screening cloud 
are, however, different from ours. We refer to some of his results below. 



A. The local spin susceptibility. 



Spin susceptibilities of the multi-channel Kondo problem also satisfy RG equation Eq.(2.21). However, the diagrams 
which contribute to the same order in 1/k are different from the single-channel case. Since the low-temperature fixed 
point for coupling constant is ^ 1/fc, each vertex produces a l//c factor. [Here we assume that the bare coupling, Ao, 
is also 0(l/k).] Each loop, on the other hand, gives a large factor of k. Combination of these factors determines the 
diagrams that one needs to calculate to a given order in l/k. The number of diagrams is finite (see Appendix A for 
details). We shall calculate the spin correlators of interest up to the first non-zero order in 1/fc. 

The solution of the scaling equations for the coupling constant up to subleading order in 1/fc were obtained by 
Gano. From the calculation of the conduction electron self-energy he finds that the /3-function is given by: 



PIXe) - -A| + ifcAl + haX% - Ifc^Al, (4.1) 

where a is some non-universal number, which depends on the cutoff procedure. The flow for the overscreened Kondo 
model is shown in Fig.|l|. The low-temperature physics is determined by the intermediate-coupling stable fixed point 
A* given by /3(A*) = 0: 



11 



2 -4a 



(4.2) 



The position of the fixed point is not universal. On the other kand, the slope of the beta-function at this fixed point, 
A = /3'(A*) is the dimension of the leading irrelevant operatoiO, and should be universal: 



A = /?'(A*) = 



k + 2' 



(4.3) 



This fact is readily checked from Eq.(|4Jj). 

It is sufficient for our purposes to consider the /?-function in the leading order in 1/fc, Eq.( |2.3| ). At this order 
A* = A = 2/fc. Solving Eq.(|2.18[), we obtain: 



^(,^£+ln|0£|) = ln^. 



where 



T, 



K 



D exp 
2 



kXf 



'k 
2 

1 



1 

a;^ 



fcAo/2 
1 - kXo/2 



k/2 



(4.4) 



(4.5) 



We assume that the bare coupling Aq is sufficiantly weak on the 1/fc scale, Aq < 2/fc. Then the solution for the 
running coupling constant is rewritten as: 



A, 



Here ^■'(y) is the function inverse to F{x) 
useful: 





' / \ A' 








+ 1 



A, 



' (4.6) 
a;exp(a;). The asymptotic form of this solution at E <^ Tx is also 



A* 



The analysis of the local spin susceptibility is parallel to the single-channel case (Section III) . It is easy to see that 
the uniform part of the local spin susceptibility should vanish in the multi-channel model as well (see Appendix B). 
Therefore, for the most general magnetic impurity (i.e. with gyromagnetic ratio not necessarily 2) we are left with 
electron and impurity parts of the oscillating local spin susceptibility X2kp,ei{'r-, T), X2kF,imp{r, T). The RG equations 
that these quantities satisfy were considered in Section II. The only difference with the single-channel case is that the 
/3-function and anomalous dimension are different, with 7imp(A) now given by: 



7iinp(A) 



fcA^ 
2 



The non- universal scale factor for the solution of the RG equations Eq.( ^.17 ) then is: 



exp 



^" 7imp(A) 

I3{X) 



dX 



1 



1 - fcAo/2' 



(4.8) 



(4.9) 



The scaling functions in the large-k hmit are determined from the perturbative analysis (see Appendix A). We find 
again that the scaling equations Eq.( 2.21 ) are obeyed, and the solutions are given by: 



X2kf.,imp{XT, ^) ^ Y 
X2kp,totiXT,x) = 



kXo/2 sinh(27ra;) 
(37r3/8)fcA|;X , X^ll(AT,a:) 



(4.10) 



sinh(2j 



sinh(27ra;) 
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where 



(^V4)(l-fcAT/2)2 



1 - fcAij/2 



(4.11) 



and \e is the coupHng at the energy scale E{ x) — T/[l — ex p(— 47 rx)]e, just hke in th e si ngle-channel case. A^;, At 
are functions of E/Tk or T/Tk given by Eq.(4.6). Using Eq.( 4.10| ) together with Eq.(4.£) we determine the scaling 
functions X'2kp,imp and X'2kp,ei up to the leading order in 1/fc in the scaling limit Ao ^ 0, r ^ l/fci? at all temperatures: 



{T/Tk, x) 
X2kF,el {T/Tk,x) 



F(-^)[(T/Tk)^] 



1 



2sinh(27ra;) \F(-^)[(T/Tk)'^] + 1 J F'^-^)[{E /Tk)^] 

Stt^x 1 
2fcsinh(27rx) (^(-i) [(£;/r^)A] j^i^' 



(4.12) 



It is interesting to note that Eq. ( 1.11 ) has a factorized form, where the T-dependence is once again that of the spin 

2Ti^Txtt{T) 



X2kp{^T,x) 



kFi-^)[{E/TK)'^] 



(4.13) 



Consider now Eq.( [1.10 ) in various limits. Obviously, A^; = Ar for r <C ^t, and A^; = At for r '3> £,t- At high 
temperatures ^t ^k, and the crossover at r ~ does not happen - just like we have seen in the single-channel 
case. For r S> ^t the correlation functions decay exponentially, just as we have seen in the single-channel case. The 
most interesting is the low-temperature limit T <C Tk, r <C ^t- In this limit we find: 



T r 

X'2kp,imp I ^ 'T~ 

T r 

X2kp,el I T~ 



(T/Tk) 



2A 



1 



4r F(-i)[(eK/[4^r])^] l-fcAo/2 
37r2 1 



(4.14) 



4fc 



T r 

ikXo/2)x2kp,im,p { 7i^,T~ 



The scaling function for the electron piece in the limit Aq — > appears in the subleading order in 1/fc. For non-zero 
bare coupling there is also a piece in the leading order, which is proportional to the impurity scaling function in 
Eq.(4.14) and the anomalous factor. 

As in the single-channel case, the weak coupling behavior is not recovered inside the screening cloud. Outside the 
screening cloud, for T ^ Tk and <C r <C Ct, the local spin susceptibility takes the form: 



X2kF,tot 



(- 


rT\ 


_ TTVF 


(AnrT^\^ 


37r2 




VF ) 


~ ArT 


KiKTl) ^ 


4fc 



(4.15) 



The T-divergence is not removed at low temperatures, and Eq.(4.15) does not have the Fermi-liquid form, as one 
could expect for a non- Fermi- liquid low-temperature fixed point. The distance dependent Knight shift for overscreened 
Kondo fixed point can also be undecstood at r > S^k using the generalization of the Nozieres' Fermi liquid approach 
developed by one of us and LudwigEZI. 

The spin susceptibility is obtained as the leading term at the low-temperature fixed point plus corrections in the 
leading irrelevant operator. For the o versc reened Kondo fixed point the leading irrelevant operator contribution 
corr espon ds to the second term in Eq.( 4.15 ). It is surprising that the dominant divergent term [the first term in 
Eq.(4.15)] in the limit T ^ 0, r ^ ^t, or in the limit fc ^ oo appears in the first order in the leading irrelevant 
coupling. An interested reader can find the details on this technical point in Appendix C. 



B. Integrated susceptibilities. 

As we have seen for the single-channel case, the static spin susceptibility is mostly given by the impurity-impurity 
correlation function, Xii{T), while other pieces contribute only a small fraction which is proportional to the bare 
coupling constant, or bare coupliag constant squared. It is easy to see that this is also the case for the multi-channel 
model. Indeed, according to Gani, 
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X^^(T) 



4T V T 



Xok 



(4.16) 



8T 



T 



Xee{T) = 



^fi-fcAgm^ 
ler V T 



Thus, from the scahng equations Eq.( ^.2l| ) we obtain: 

Xtt 



Xr. 



Xtt = 



Xtt 



(l-fcAo/2)2 l-kXo/2' 
where the scaling function for the total spin susceptibility is given by: 



Xtt 



kXn 



Using Eq.(4.6), we can rewrite Eq.(4.18) in the form: 



1 / i^(-i)[(T/r;^)^] 



Xtt{T) \^^(-i)[(y/y^)A] + I 

The electron-impurity and impurity-impurity correlators contain smallness associated with the bare coupling: 



Xi' 



kXo , s 

'—Xtt[T) Xe 



ex, 



°-Xtt{T). 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



Thus, the spin susceptibility is given mostly by the impurity-impurity spin correlator, and for a system with impurity 
g-factor 5^2 there are corrections to the bulk susceptibility proportional to the bare coupling. For T < Tk the 
scaling function for the total spin susceptibility takes the form: 



1 



Xtt ~ , 



T 



2A 



(4.21) 



As in Section IIIB, this fact is easily understood in the bosonic language using canonical transformation. The bosonized 
Kondo Hamiltonian for the fc-channel model has the form: 



H 



dr 



heVk 



dripL{r) 



+ i/P"™ + V2^XoS'dr^L{Q) 



+const ■ Ao(5+e*v8V^'^^.(o)c)?"^™ + h.c.) - h,S' 



(4.22) 



Here ip is the canonically normalized total spin boson, i.e. the sum of the spin bosons for each channel divided by 
^/k. The additional, independent degrees of freedom which couple to the impurity correspiMid to the SU{2)k Wess- 
Zumino-Witten model with one free boson factored out. This is the Zk parafermion modelEZI. For the k = 2 case it 
corresponds to an extra Ising degree of freedom, or equivalently a Majorana fermion. These extra degrees of freedom 
play no role in the canonical transformation. 

Changing the bosonic field pl{i') — <^l(^) + '^j^'^ , the Hamiltonian takes the form: 



H = 



k 



dr{dr^L{r)f + -ffo"'"'" - ^hl{2L) + \/2^Ao5^9,(^l(0) 



+const ■ Ao(S•+e*^/8^'^^(o)OP"™ + h.c.) - (^,kXa^ 
Thus, for the Free energy, we have: 

hekXo-. 



T k 

F{K,h,)^-- hl + F{0,h 



-) 



(4.23) 



(4.24) 
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For the susceptibilities we then obtain: 



- - f (A 

Xie — 2 Xiii Xee — \ 2 / \ 

with xo = 27™7^' Pauh term. This agrees with the large-k results. 

Let's now return to the issue of screening. The electron-total piece of the spin susceptibility, Xet = Xee + Xei, is 
given by the integral of the local spin susceptibility xC^iT")- As in the single-channel case, since x(r, T) only has the 
oscillating piece, this integral is determined by the short-distance contribution, r ^ l/kp- The form of x(^j'^) 
r ~ 1/kp is cutoff-dependent. However this dependence disappears in the integral, which describes conduction electron 
spin polarization. In case of a 3D Fermi gas the cutoff procedure is well-defined. The fact that the net conduction 
electron spin polarization due t o imp urity comes mainly from r ~ l/kp is indeed justified to the orders we worked in 



perturbation theory. From Eq.(4.10), with Xe — — Ao ^ 1/k we can write for the local spin susceptibility: 



We have checked this conjecture to the leading order in 1/k. Integration of this expression over r gives the correct 
result for Xet, 

Obviously, the major contribution to the integral 

/■°° „ , , f°° d /sin2fci.r\ 

comes from r ~ l/kp- 



C. Equal-time correlation function. 



As we have seen above, the zero-frequency spin correlator vanishes as T^'^ when T ^ 0. It also obeys a zero-sum 
rule. As in the single-channel case, the equal-time spin correlator K{r,T) — (^^((r, O)S'f^p(O)) has a nonvanishing 

sum rule, since (S^iSfj^p) = —1/4. The uniform part of the equal-time spin correlator is non-zero. 

Consider th e eq ual-time correlators i^„„(r, T) and K2kF{f',T) using the \/k expansion. K(r,T) satisfy scahng 
equations Eq.( 2.21 ). As in the single-channel case, for r > the spin correlators decay exponentially. The behavior 
is most interesting for r <C Ct, where our expressions are considerably simplified. Expressing our results in terms of 
effective coupling at scale r, A^, we get: 



Kun{K,r/^T) 

K2kp{K,r/S,T) 



1 TTVpkXl 

1 - fcAo/2 2^ 

1 TTVp 



(1 - /jA^/2)2 



(4.29) 



l-fcAo/2 8r 



kXr{l - kXr/2) 



We can rewrite these expressions using Eq.(4^) in terms of T/Tk, t/S^t variables. Suppressing the anomalous factor 
1/(1 — kXo/2), we obtain: 



T r 

K2kF ( ^ 



kr 



(4.30) 



4r 



L[{ATir/UY 



where L{x) is the function defined by 



L{x) 



F(-i)(l/x) 
(F(-i)(l/x) + l)2 



(4.31) 
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A plot of this function is s hown in Fig. ^. As we have discussed in Section III, the integral of K{r, T) should not 
vanish. It is given by Eq. (|3.35 ), as in the single-channel case. The integ ral o ver long distances r ~ can be 
calculated exphcitly from Eq.(4.2£) by changing variable r ^ Xr- Using Eq.(2.18), 



drKun{K,r) 
2'kvf 



A* 



2-kXo2X* 



kXo 



(4.32) 



Thus, in this case the screening length ^ S,k is explicitly present. The de pend ence on the bare coupling constant 
Ao is surprising, since it should not be there according to the sum rule Eq.( 3.35| ). The missing part of the sum rule 
comes from the short distances. To provide the most transparent demonstration of this, we write the seco nd eq uation 
in Eq.(4.29) for a 3D Fermi gas, so that cos{2kpr) is replaced by cos(2/cpr) — [sm(2kpr)/2kpr], as in Eq.(4.26). The 
short-distanc e inte gral, which is analogous to Eq.( 4.28| ), gives precisely the compensating term — /cAq/S needed for the 
sum rule Eq.( 3.35 ) to be obeyed. 

The low-temperature decay of the equal-time correlator at r ^ in the overscreened multichannel Kondo model 
can be obtained using conformal field theory approach (see Appendix C). Indeed, at the low temperature fixed point 
we haveO'El: 



const(f){0,0)Tj^'^, 



(4.33) 



where cf) is the s = 1 primary of dimension A = const is a non-universal constant. We then obtain for K; 



from conformal invariance: 



2kf 



K,,Ar) cx (^^^^^V.1(0,.)^.(0,-.)) « 



(4.34) 



in agreement with the large-k result Eq.( 4.30| ). The same leading order calculation gives zero for Kum since 
(^(0, 0) • J(0, r)) = 0. The first-order term in the leading irrelevant operator gives 



Kun{r) oc 



2T2Ay,l+2A ' 



(4.35) 



which also agrees with Eq.( 4.30| ). It is interesting to note that, unlike the single-channel Kondo model, the long- 
distance decay of the uniform and 2kp correlators is different. 



V. CONCLUSION 



Although the techniques employed in this paper, renormalization group improved perturbation theory and the large 
k limit are of limited validity, they have led to one exact result (all orders in perturbation theory) and suggested a 
certain conjecture which, if true, lead to a rather complete picture of the Kondo screening cloud. We first summarize 
the exact result and the conjecture, pointing out a consistency check between them and then state the resulting 
conclusions. 

• The uniform part of the r-dependent susceptibility, vanishes to all orders in perturbation theory. On the other 
hand, the equal time correlation function has a non-zero uniform part, varying on the scale at T = 0. 

• The 2kp part of the r-dependent susceptibility has a factorized form at vp/T ^ r: 

X2kAr,T)^ f{r)xu{T), (5.1) 



where Xtt{T) is the total susceptibility (less the free electron Pauli part). At small r this becomes: 

X2kAr,T)~^ ^XrXttiT), (5.2) 



where is the effective coupling at scale r. This was verified to third order in perturbation theory and in the 
large k limit (including the 0{l/k) correction). 



16 



There is an important consistency check relating this result and conjecture and the result Xie = ^(-^o/2)xiij 
following from the formula: 

Xet = j d\x{r). (5.3) 

Since xC*") is an RG invariant, it has no explicit dependence on the bare coupling. If the uniform part had been 
non-zero, its integral would have given a contribution to Xet which would be unsuppressed by any powers of the bare 
coupling. The integral involving X2fcF(^) gives for r ^ 1/fcF due to the cos{2kpr) factor and hence is determined 



by the value of X2kp at short distances ofO(l/fcF). In this limit X2fcF(^) oc Ar ~ Aq and integrating Eq. (5_^) gives 
Xet ^ -(Aofc/2)xtt. 

Strictly speaking this consistency check requires yet another conjecture: 



X{r) 



Vp 



, sm2kpr 
cos[2kpr) 



2kpr 



kXrXttjT) 
87rr2 



, sm2kpr 
'^"^(2^^'^) - -2k^ 



(5.4) 



for r <C ^K,vp/T. This last conjecture, involves corrections of 0(1/ kp) which we have not calculated systematically 
and go beyond the scope of the one-dimensional model. We did check the result in lowest order in l/k. 

Despite the limitations of our calculational approach, we are thus led to a fairly complete understanding of the 
Kondo screening cloud. The heuristic picture of Nozieres and others of the Kondo groundstate is seen to be correct. 
The impurity essentially forms a singlet with an electron which is in a wave-function spread out over a distance of 
0{^k)- This is seen from our calculation of the T = equal time correlation function which varies over the scale ^k- 

On the other hand the behaviour of static susceptibilities is considerably more subtle. A naive picture that an 
infinitesimal magnetic field fully polarizes the impurity but induces a compensating polarization of the electrons is 
certainly wrong. Rather the impurity polarization is proportional to the weak magnetic field and the integrated 
polarization of the electrons (with the free electron value subtracted) is much smaller (proportional to Aq). The 
finiteness of the T — impurity susceptibility results from its tendency to form a singlet with the electrons. 

If we now examine the r dependence of the electron polarization, we find that it is small at short distances (0(Ao)). 
However, it exhibits a universal oscillating form at long distances which is not suppressed by any powers of Aq but 
only by a dimensional factors of l/r^. The fact that it is purely oscillating ensures that the contribution to the 
integrated polarization is negligible. The envelope of this oscillating susceptibility, consists of the dimensional factor 
of l/r^ times an interesting and universal scaling function of r/£,K and T/Tk- This scaling function factorizes into 
Xtt{T)f{r/^K) for vp/T:^r. 

Our work leaves various open questions for further study. It seems plausible that our conjecture could be proven 
to all orders in perturbation theory, thus putting this work on a more solid foundation. There are three interesting 
universal scaling functions which we have introduced, one for the 2kp susceptibility and two for the uniform and 2kp 
equal time correlation functions. A general calculation of these functions could perhaps be accomplished by quantum 
Monte Carlo or exact integrability methods. Results on the T ~ limit of the susceptibility scaling function were 
given in Ref. [|^. An obvious generalization of our calculations is to general frequency dependent Green's functions. 

Most importantly, experimental results on the Kondo screening cloud are very limited. The NMR experiments of 
Boyce and Slichter only probe extremely short distances, r « 1/kp. Our work shows that these results are entirely 
consistent with a large screening cloud. However, these experiments do not directly probe the scale £,k- NMR is 
probably not a feasible technique for doing this since it is difficult to study distances of more than a few lattice 
constants. One possibility might be neutron scattering, which_could in principle measure x{Qt^) for Q ~ 2kp. An 
alternative is to study small samples with dimensions of 0(^if )E3. 

We would like to thank A. V. Balatsky, J. Gan, F. Lesage, N. Prokof'ev, H. Saleur, D. J. Scalapino, E. S. S0rensen, 
P. C. E. Stamp, B. Stojkovic, C. Varma and A. Zawadowski for useful discussions and comments. This research was 
supported by NSERC of Canada. 



APPENDIX A: PERTURBATIVE RESULTS 



The diagram technique for interactions involving spin operators is complicated due toJjheir nontrivial commutation 
relations. It is possible to express these operators in terms of pseudofermion operatorsE3cj: 

S™p = ^ E f^'^'^afp- (Al) 

a, 13=1,2 
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The problem in using the fermion substitution Eq.(Al) is that the cr-matrices have dimensionahty 2, while the fermion 
space is four-dimensional. Thus, only the states with 



N = Y.f^"f- = l (A2) 

a 

are physical. This constraint is imposed by choosing appropriate chemical potential^. For example, Popov's 
techniquecj adds an imaginary chemical potential, z7rT/2, to the pseudofermions. Then the contribution of the 
nonphysical states to the partition function is zero. The diagram technique then becomes the standard fermion 
technique with the one-dimensional conduction electron (left-movers) propagator {iujn + VFk)~^ , the pseudofermion 
propagator {iuJn ~ , and the interaction Hamiltonian 

H^nt = VFXo^^^^i"iO)p-'f5^L0iO). (A3) 

For our purpose of computation of spatial correlators it is convenient to work in the coordinate r, r- space, where the 
propagator for the left-movers takes the form: 

ttT 

Goiz) ^ . z^vpT + ix. (A4) 

sm [ttI z\ 

For the lowest-order diagrams it may be more convenient to calculate time-ordered impurity spin averages directly. 
Such spin operator Green's function approach was applied successfully, for example, in case of long-range Heisenberg 
ferromagnetSLj. Consider 

{S'^S^ ....S^) , (A5) 

where i,j,...,k — }. Obviously, this average is zero when the total number of operators is not equal to 

the total number of . Consider first averages containing only 5*^ operators. For odd number of spin operators it 
vanishes. In our simple S* = 1/2 case Tr ([5^]^") = 1/4". One can use spin commutation relations and the relations 
S+S- = (1/2) + 5^ S-S+ = (1/2) - 5^ to calculate the average Eq.(^). 

All diagrams for the spin susceptibility x(r, T) up to third order are shown in FigJ^. The graphs (a)-(d) represent the 
electron-impurity part, while the graphs (e)-(i) the electron-electron part. We only show the electron Green functions 
on these diagrams. The dashed line represents the boundary. For the electron-impurity spin correlation function the 
external electron spin operator Sei (r) takes the propogator away from the boundary. In case of electron-electron part 
of the Knight shift there are two such operators. We have to integrate over t he po sition of one of these operators. 

Straightforward calculations lead to the final results stated in Eqs. (3.1),(^T7) of Section IIIA. To calculate the 



equal-time correlator {Sei{r, 0)Simp{0)), we need to evaluate the graphs (b)-(d) of Fig. once again. The first graph 
(a) is frequency-independent, i.e. it is the same as for the electron-impurity part of the local spin susceptibility. Both 
uniform and 2kF parts are now non-zero. The result of this calculation is given by Eq.(|3!36|) of Section IIIC. 



For the discussion of static susceptibilities in Section IIIB wc need to calculate impurity- impurity part, in addition 
to space integrals of Xieii^jT) and Xeeii^jT). The second- and third-order graphs for Xii{T) are shown in Fig.|^. The 
leading order is, of course, 1/4T. We find that 

ATxu^l~Xl(\n^+A,^ (A6) 

- Xl f In^ — + ^2 In — + const 
ATx,e = - Y + SiAg + ^ In ^ + constXl 
4rXee = -^+const\l. 

In general, the constants Ai, A2, and Bi in Eq. ([A^) depend on the cutoff procedrpe. However, these three constants 
are connected, A2 + iBi — 2Ai = 0, as follows directly from the results of Wilsoro on the scaling properties of the 



total spin susceptibility. Using this connection and Eq.(3.25), the fact that all three scaling functions for the spin 



susceptibility are equal up to the terms ^ Aj^ is easily demonstrated. 

Consider now the multi-channel case. As we have mentioned in the text, the graph selection in this case is different, 
since each vertex A is ^ 1/fc. To the order 1/fc we need to calculate all the graphs in Fig. |^, except (c) and (i), which 
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are of the order In additi on, we need to calculate the fourth-order graph shown in Fig. ^. The result of this 

calculation is given by Eg. ( 4. 10 ) in the text. 

Calculations of the equal-time correlator are somewhat more involved. While X2feF('', in Eq.(4.29) is also non- 
zero up to this order, Kun{T, T) vanishes. We need to go to the next order in 1/fc to find the answer. For the terms of 



the order 1/fc^, we need to calculate graph (c) in Fig 

The bulk susceptibility results are found again by calculating x 
leading order we only need to consider second-order graph in Fig. p 



and additional fourth and fifth order graphs shown in Fig |10| 
(T) and r- integrating the Knight shift. In the 



APPENDIX B: PROOF THAT THE UNIFORM PART OF THE LOCAL SUSCEPTIBILITY VANISHES 



A s cla rified in the text, the local spin susceptibility can be written as a sum of impurity and electron parts (see 
Eq.( ^.lO )). We will consider these two parts separately fo r the purpose of this proof. 
Consider first the impurity part, Xun,imp(''')- Using Eq.(2.15), one can write: 



Xun,h-npir,T) = VF 

+ (r ^ — r) 



(Bl) 



where Hint is given by Eq.( |2.9| ). The fact that this contribution vanishes is very easily seen when we perform the r 
integration. Indeed, in every order in perturbation theory Xun,irap{i',T) can be written as: 



/ / dTidT2l{Ti,T2,r) y. F{ti,T2) 
Jo Jo 



where 



or, equivalently, 



Jo 



)G{-r,T2-T), 



1 = 



dT{TlTf 



Q sm['KT{vFT — vpTi + ir)] sm[T:T{vFT2 — vft — ir)\ 



(B2) 



(B3) 



(B4) 



After the change of integration variable, r — s- ex'piilixTvF'r), one encounters contour integration with two poles on 
one side (see Figjll|) , and / = 0. 

Consider now the electron part, Xttn,ci('')- Here the cancellation of x«n.ci('') is less trivial since there are other 
graphs in addition to those with the integration Eq.(p^ (see Fig.p^: 

G(r', -Ti)G(r - r', T)G(-r, T2 - t)t]{t2 - n) (B5) 
G(r, T - T\)G{r' - r, -r)G(-r', T2)?7(t2 - n), 

where r\{T2 — t\) is determined by the full perturbative series. We now introduce the complex notation, z = TrT{vFT + 
ir), and remember that G(z) = ttT/ sinz. Then the sum of the graphs in Fig.|l^ gives: 



V{z2 - zi) 
sin(z — z') 



sin(z2 — z) sin(z' — zi) sin(z2 — z') sin(z — zi) 

'n{z2 - zi)sin(zi - Z2) 

sin(z — zi) sin(z2 — z) sin(z' — zi) sin(z2 — z') ' 



(B6) 



which is graphically presented in Fig.|l2|. Integration over r Eq.(B4) yields zero in this case as well. Generalization of 
this proof to the multiple number of channels is quite trivial. Indeed, the graphs that cancel have the same channel- 
dependent factor. As we have seen above, the crucial step of the proof is that the integral Eq.(B4) is zero. Thus, the 
g = part of the correlator is absent only for zero-frequency spin-spin correlators, not for equal time correlators. 

Note that the absence of the uniform part in the distance-dependent Knight shift becomes trivial in the bosonic 
language (see Section IIIB). Indeed, since 
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SliivpT + ir) = —=dr(pLivFT + ir), (B7) 
V 27r 



we find: 

Xun{r) ^ [ dr l^dr(pL{vFT + ir)S'A = ^(< (pL{ir + vfP)S' > - < (pL{ir)S'' >) = 0, (B8) 

Jo \v27r / v27r 

because ipL^z) is periodic in the imaginary time variable. Note that we don't need to worry about a potential short- 
distance singularity because the total spin has been replaced by the impurity spin in the expression for Xun using the 
above argument. A similar argument for Xun = was given in Ref . [ [l3| . 



APPENDIX C: LOW-TEMPERATURE LONG DISTANCE LOCAL SUSCEPTIBILITY IN THE 

MULTI-CHANNEL KONDO MODEL. 

X2kF{r,T) is determined by the infrared stable fixed point for r 3> Ck, T <^ Tk and any value of the ratio rT/vp. 
For k > 1 (and Simp = 1/2) this fixed point is of non- Fermi liquid ty[pej-The low-temperature non- Fermi- liquid multi- 
channel Kondo fixed point was analysed _bv Ludwig and one of ual3lla using conformal field theory. We refer the 
reader to these works and a recent reviewEll for details. In the bosonized form spin, charge, and flavour sectors of the 
free fermion Hamiltonian are separate. Only the spin sector is interesting in the-J^ondo problem, since the impurity 
spin couples to the spin current. The effect of the strong coupling fixed pointEJ is such that the low-temperature 
Hamiltonian density is written in terms of new spin currents. 



where 



3{x) ^Y.^{^{x)^^Lj{x) + 27rSS{x). (C2) 



The Fourier modes of the spin currents for a system with Hamiltonian density Eq.(|Cl|) defined on a large circle of 
circumference 21, 



J„ = / da;e™^-/'j(a;), (C3) 
27r ./„/ 



satisfy the usual Kac-Moody commutation relations, 

[j;^, Ji] = le^'^J^l+m + ^knS^''5„+m,o- (C4) 

Here e'^'"^ is the antisymmetric tensor and k is the Kac-Moody level. To the leading order, the Knight shift is given 
by: 

X2k, (r, T)^-^J^^ drdy (^^[(0, r)^^i(0, -r) J^(r, y)^ . (C5) 

Using OPE 

^iQ^^L{z) = -ll^^^{z) + Regie (C6) 
2 C- z 

^V'l(^) = ^^I(^) + R^9iC - z), (C7) 
where Reg{C, — z) denotes a function which is regular at ^ z, we rewrite X2kF (^j T = 0) as 

<VL(0,-r)Vl(0,r)>. (C8) 



X2kj,{r) = / / drdy 



oo <J —oo 



T + iy — ir t + iy + ir 
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where 



The Green's function for two points on the opposite sides of the boundary takes the form: 

^l(zi)^i(^-2)) - (C9) 
/ 2:1-^2 

cos[27r/(2 + fc)] 

- cos[V(2 + fc)] ^^^"^ 

is the S- (scattering) matrix, calculated in Ref.[ This is a universal complex number, which depends on the 

universality class of the boundary conditions. In the one-channel Kondo effect S'(i) = — 1, corresponding to a 7r/2 
phase shift. At the overscreened Kondo fixed points |5'(i)| < 1, which means multiparticle scattering. Substracting 
free electron contribution and performing the integrals, we find: 

X2kAr) = k^—^- (Cll) 
In the limit fc — s- 00 this gives X'2kp{'>') — 37r^ /4fc, in agreement with the large-k result of Section IV. For fc — 1 



it agrees with the Fermi liquid result Eq.(3.15). Note that no anomalous power laws occur in the leading order in 
irrelevant coupling constants. Only the normalization reflects the non-Fermi liquid behavior. As in the single-channel 
case, finite temperature calculations multiply this expression by the factor 27ra;/ sin(27ra;), where rT/vp- 

Consider now corrections to this expression. The leading irrelevant operator which appearsliZI in the effective 
lagrangian at the overscreened Kondo fixed point is J_i • cj), where </> is the s = 1 SU{2) KM primary field with the 
dimension A = 2/(2 + fc). The dimension of this singlet operator is 1 -I- A. We can again write this additional piece as 

H,nt-T^iJ-l-m)- (C12) 



Thus the correction is given by: 



5x2kAr,T) = -,^ / / drdridy^V^i(0,r)— ^i(O,-r)J-(T,2/)(J_i.0(ri,O))). (C13) 



VF 

To find the most singular part of this expression as r ^ 0, we use the boundary OPE 



2 

From conformal invariance, this zero-temperature correlator 



^i{0^^r)^MO-^r)^^^. (C14) 



((0(0) • J(zi))(J-i • 0(Z2))) = I 1,^1 ^' 1,^ . (C15) 



The finite-temperature correlation function which appears under the integral in Eq.(C13) can be obtained using 
conformal mapping, a conformal transformation which maps the finite-temperature geometry (half-cylinder) onto the 
zero-temperature half-plane. 

z = tan(7rTw). (C16) 

Here w ^ t + ir ui the finite-temperature geometry. A Virasoro primary operator A{z) of left scaling dimension A^i 
transforms as 

^H=(^) A{z), (C17) 

under conformal transformation. Using dw(z)/dz = l/7rr(l + z^), we express the finite-temperature correlators in 
terms of the zero-temperature ones. The net effect is such that the fact ors 1 /(2:1 — Z2) for the half-plane get replaced 



by ttT/ sin(7rT[i(;i — W2\) on the half- cylinder. Doing the integral in Eg. ( CIS ) and dropping the constants, we obtain: 
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in agreemen t wit h the large-k result of Section IV, This term is subdominant, for r <C vp /T, compared to the leading 



term in Eq.(Cll). On the other hand, it becomes larger than the "leading" term if we take T ^ with r ^ S^k held 



fixed. Anomalous powers appear from irrelevant operator corrections. 
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FIG. 1. RG flows for the single-channel and the multi-channel Kondo problems. 
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FIG. 2. Singular third-order graph for x{f,T). 
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FIG. 3. Scaling regimes for x^kpi^TjX = r/^r) in the single-channel Kondo effect. 
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FIG. 4. Scaling regimes for the oscillating part of the equal-time spin-spin correlator K2kp (At, x = r/^r) in the single-channel 
Kondo effect. 
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FIG. 5. Scaling regimes for the uniform part of the equal-time correlator if„„(Ar,a; = r/^r) in the single-channel Kondo 
effect. 
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FIG. 11. Contour of integration for I. 
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FIG. 12. Cancellation of the uniform part of the local spin susceptibility. 
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